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New BABAR measurements are presented for the exclusive rare decays B — > K^*'t^f.~ including branching 
fractions, isospin asymmetries, direct CP violation, and lepton flavor universality for dilepton masses below and 
above the J/i/i resonance. Unexpectedly large isospin asymmetries are observed in both K£'^£~ and K''£^£~ 
decays. For the combined K£'^l~ and K*£~^£~ data a 3.9cr significant deviation from the SM prediction is 
found. Furthermore, recent BABAR results from an angular analysis in B — > K*t^l~ are reported in which 
both the K* longitudinal polarization and the lepton forward-backward asymmetry are measured for dilepton 
masses below and above the J/i/i resonance. Finally, results of recent searches for B — > 7ri'+^~ from Belle and 
B — > K^*'i'P from BABAR are summarized. 



1. Introduction 

The flavor-changing neutral-current (FCNC) pro- 
cesses, b —^ s{d)i£, provide an interesting hunting 
ground to look for new-physics phenomena. In the 
Standard Model (SM), FCNC transitions are forbid- 
den at tree level, but they are allowed to proceed 
via electroweak-loop and weak-box diagrams. Lowest- 
order processes are depicted in Figure [TJ The elec- 
troweak loops consist of a W-boson and at, c or 
M-quark, where the i-quark contribution dominates. 
One of the loop particles radiates a photon or Z-boson 
to conserve momentum. In the weak-box diagram the 
b to s{d) transition occurs via the emission of two 
W bosons. The emitted lepton pair is e^e~, ^'^ fi^ , 
T~^T~, or I'D. Final states with a t^t^ pair, how- 
ever, are difficult to measure because of the missing 
neutrinos and are not discussed further. 

In chapter 2, we discuss properties of the exclu- 
sive decays B — > K'^*H^l~ where t^(.~ is either e+e^ 
or fJ-'^ n~ and introduce the observables we measure. 
In chapter 3, we focus on BaBAR measurements of 
B -^ K^*'£~^i~ branching fractions, isospin asymme- 
tries, direct CP violation and lepton flavor universality 
after discussing the analysis strategy. In chapter 4, we 
discuss the BaBAR angular analysis of B ^ K*i^£~ 
decays. In chapter 5, we present the latest branching 
fraction upper limits oi B —^ TT£~^i~ decays. In chap- 
ter 6, we show results of recent searches for B -^ Ki/D 
and B —* K*^^ before concluding in chapter 7. 



2. Properties of 5 ^ K(*)£+£ Decays 

The SM calculations are based on a low-energy ef- 
fective Hamiltonian [l|] that factorizes perturbatively- 
calculable short-distance contributions expressed by 
Wilson coefficients, Ci, from the non-perturbative 
long-distance contributions of the b ^ s transition 
operators, Oi. In lowest order, dominant contribu- 
tions result from the magnetic dipole operator Oy, 
the vector-current operator Og and the axial-vector- 
current operator Oio- While Oy represents the pho- 
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Figure 1: Lowest-order SM processes for b -^ s{d)££; (left) 
electroweak penguin loops and (right) weak box diagram. 
If £ is a charged lepton (neutrino), £' is a neutrino (charged 
lepton). 



ton penguin diagram, Og and Oio result from linear 
combinations of the weak penguin and weak box dia- 
grams. Thus, the relevant Wilson coefficients for these 
modes are C7, Cg and Ciq. QCD effects, however, 
introduce operator mixing. Besides as corrections, 
the Wilson coefficients also receive contributions from 
other operators [3, l3|, y|. It is customary to absorb 
these by defining effective Wilson coefficients Ci. The 
Wilson coefficients are functions of the renormaliza- 
tion scale /j, and the squared dilepton mass, q^. For 
low values of q^, /i = 4.6 GeV and a top quark mass 
fht{mt) — 167 GeV/c^ in the MS renormalization 
scheme, the SM predictions at next-to-next-to-leading 
order (NNLO) yield Cy = -0.31, Cg =jl.21 and 
Cio = -4.31 d S, [1. The magnitude of C7 is con- 
strained experimentaUy by B{B — > ^^7) [6]. 

Since FCNC transitions are suppressed in the SM, 
loops and box diagrams from processes beyond the 
SM may yield non-negligible contributions. For exam- 
ple, new contributions may arise from loops contain- 
ing a charged Higgs boson or supersymmetric (SUSY) 
particles [7|, [^, [^ . Examples are shown in Figure [2] 
New physics contributions will modify the effective 
Wilson coefficients from their SM expectations. In 
addition, scalar and pseudoscalar couplings that are 
absent in the SM may modify the decay rate [j|. In 
order to have the ability for uncovering new physics 
phenomena, the SM predictions need to have suffi- 
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cient precision. Most recent calculations use QCD fac- 
torization to separate the short-distance physics from 
the long-distance effects [a, [T^]- The effective Wil- 
son coefficients entering the short-distance part are 
calculated in the NNLO approximation. The long- 
distance effects including the hadronization process 
are expressed in terms of hadronic matrix elements 
of the & — *■ s transition operators between the initial 
B and the K''*'^ final states. Since the hadronic ma- 
trix elements cannot be calculated from first princi- 
ples, they are parameterized in terms of form factors 
[llj that are calculated with the help of light-cone 
sum rules (LCSR) [12i | or in soft-collinear effective 
theory (SCET) [l|]. Though the form-factor calcula- 
tions include next-to-leading order (NLO) QCD cor- 
rections they bear large theoretical uncertainties that 
are presently the dominant uncertainties in the SM 
predictions for exclusive decays. Thus, it is important 
to measure many different observables in different in- 
clusive and exclusive electroweak penguin processes 
in order to extract meaningful results for moduli and 
arguments of the effective Wilson coefficients. 
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Figure 3: Predicted B — > K*t^l~ branching fraction in 
the SM as a function of (f with (dark curve) and without 
(light curve) contributions from charmonium resonances 
[ill ]. The shaded region indicates the SM uncertainties. 
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Figure 2: Examples of new-physics contributions in loop 
processes, (top left) a charged Higgs boson with u, c, t 
quarks, (top right) a chargino with u, c, t squarks, and 
(bottom) a neutralino or gaugino with d, s, b squarks. 

Figure [3] shows the q^ dependence of the B -^ 
K*e+e- branching fraction in the SM [ui- The^over- 
all shape is determined by the q^ dependence of C9(g^) 
except for the low q^ region that is donrinated by the 
l/q"^ term originating from B — > K*j. The singularity 
at g^ — 0, which is not present in B ^ K£^£~ , is cut 
off by the finite mg masses. Thus, B{B -^ K*e'^e~) 
is expected to be ^ 25% larger in the SM than 
B{B —> K* fjL^ n~), where the increase in branching 
fraction comes just from the extended q^ < 4to^ re- 
gion. In addition, the hadronic decays B -^ J/tpK* 
and B -^ ^{2S)K* with J/V', ?/'(25) -^ e+£~ interfere 
with signal modes. Since the branching fractions of 
the charmonium modes are more than two orders of 
magnitude larger than those of signal modes, we need 
to exclude sufficiently large q^ regions in the vicinity 
of the charmonium resonances to remove them (see 
section [3]) . 

Besides partial and total branching fractions, asym- 
metries are of great interest in testing the SM, in par- 



ticular for exclusive decays since many uncertainties in 
both predictions and measurements cancel 14]. For 
B -^ Ki~^g~ and B -^ K*£'^£~ we have performed 
new measurements of the isospin asymmetries, direct 
CP violation and lepton flavor universality. 

The CP-averaged isospin asymmetry is defined as 



Af 



B{B^ -^ K^*)^£+£-) + r6(B± -^ K(*)^£+£~) 



(1) 



where r = ro/T+ = 1/(1.071 ± 0.009) ^ is the ratio 
of 5° and B'^ lifetimes. Figure 2] shows the q^ depen- 
dence of Af' in the SM. For g^ — + the isospin asym- 
metry is 6 — 13%. With increasing q^, Aj decreases 
approaching a value ~ —1% for q^ > 2.5 GeV^/c* @]. 
The direct CP asymmetries 
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may provide useful constraints on non-SM 
physics [la] ■ The SM Acp predictions for modes stud- 
ied here are very small, 0(10"'^), and new physics 
at the electroweak scale may provide significant 
enhancements |17| . 

The ratios of rates to /i^/i^ and e+e^ final states 
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Figure 4: Predicted q^ dependence of isospin asymmetry 
Ai{q^) in the SM in the low q^ region [3]. 



are sensitive to the presence of a neutral SUSY Higgs 
boson 0]. In the SM, TZk is expected to be unity 
modulo a small correction accounting for differences 
in phase space [1|. TZk* should also be close to 
unity for mu > 2m ^. Due to the 1/g^ dependence 
of the photon penguin contribution, however, there 
is a significant rate enhancement in K*e^e^ modes 
for mg+g- < 2m^. The expected SM value of TZk* 
including this region is 0.75. In order to test this pre- 
dicted rate enhancement, we fit the K* dataset in the 
entire q^ and in the low q^ region with and without 
inclusion of events in the q^ < 0.1 GeV^/c'' region. 
Theoretical uncertainties for TZki*) predictions in the 
SM are just a few percent. For example, in two-Higgs- 
doublet models the presence of a SUSY Higgs might 
give ~ 10% corrections to TZki.*) for large tan/3 [^. 

Besides q^ , the B -^ K*t^t'' differential decay rate 
also depends on three angles, Ok^Oi, and X: where Ok 
is the angle between the K and the B momenta in the 
K* rest frame, 9i is the angle between the £^ {£^ ) and 
the B(B) momentum in the i^i~ rest frame, and x is 
the angle between the K* and £~^i~ decay planes. The 
present data sample is not large enough to perform a 
full three-dimensional angular analysis. Thus, we fit 
the one-dimensional angular distributions [IQ] 



W{cos9k) - -TLCOs''9K + -{l-TL)smHK, 
4 o 



+ Afb cos 



(4) 



The parameter J-^ is the K* longitudinal polar- 
ization. The second parameter, Afb, is the lepton 
forward-backward asymmetry. Both parameters are 
functions of q^ . Figures [H] shows the q^ dependence of 
the K* longitudinal polarization. In addition to the 
SM prediction, the distribution for the flipped-sign 6*7 
model is shown for which the sign of C7 is opposite to 
that in the SM. 
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Figure 5: Predicted q^ dependence of the K* longitudi- 
nal polarization in the SM (solid) and the flipped-sign C7 
model (dashed). 



The q^ dependence of lepton forward-backward 
asymmetry is given by 



d^ 



FB 



dq^ 



ex ^ne[{Cr>{q'))Cio{q')V(q')A,{q') 



+ ^^Cr{q^)CMq')iV{q')T2iq') 



(1 



Mb 



Ai(g2)Ti(q2)(l 



Mb 



■))]• (5) 



The functions Ai{q'^),V{q'^),Ti{q'^) and T2{q'^) are 
g^-dependent form factors. The shape of the lepton 
forward-backward asymmetry results from different q^ 
dependences of the interference terms — C9(q^)Cio((7^) 
and ~C7{q^)Cw{q^) / q^ . In the SM, the first (second) 
term is positive (negative) and increases (decreases) 
with q^ . Thus, Afb is negative at low q^ , crosses zero 
and remains positive for large values of q^ [ll|, [1^, [3 ■ 
For B -^ K*t^£~ the zero crossing is predicted in 
NLO at ql = 4.2 ± 0.6. Figure [6] shows the low q^ 
region for B -^ K*£~^i~ ^. The uncertainty of q^ is 
dominated by uncertainties in the form-factor calcu- 
lations. This is different for inclusive modes. Recent 
calculations of Afb in i? ^ Xsi~^£~ at order 0(0^) 
(NNLO) that include in addition electromagnetic cor- 
rections predict a g^ dependence shown in Figure [7] 
|19| . The zero crossing is expected at ^q = 3.5±0.12. 

Figure [5] shows the SM prediction of Afb{q^) in 
the entire q^ region. New physics may change both 
magnitude and phase of the Wilson coefficients [23| 
yielding shapes that differ from those in the SM. 
We consider three simple examples, also shown in 
Figure [51 in which the signs of the two interfer- 
ence terms are reversed (i.e. a phase change by tt) 
P Q m, HH [13, m, m. changing the sign in C7, 
yields a positive Afb for all values of q^, whereas a 
reversed sign in C9{q^)Cio{q^) yields a negative Afb 
in the entire q"^ region. If the sign of both interfer- 
ence terms is flipped, the resulting distribution has 
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Figure 6: The low q^ region of Afb for B -^ K'l'^l 
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Figure 7: The low g^ region of Afb for B -^ Xsi+T [H 



the negative mirror image of the SM distribution. 



3. Study of S ^ K'-*H+i~ 

We have studied the exclusive decays B -^ Ki~^£~ 
and B -^ K*i+£~ using 384 milhon BB pairs col- 
lected at the T{4:S) resonance with the BaBAR detec- 
tor [23] at the PEP-II asymmetric-energy e+e" col- 
lider. We reconstruct ten individual final states, in 
which a K+, K0(-> tt+tt"), K*°(892)(-> K+tt") or a 
K*+(892)(-> K+7r° or K°7r+) recoils against an e+e" 
or a /^^/i~ pair. We require good particle identifica- 
tion (PID) for e^,iJ,^,K^ and tt^. We select elec- 
trons (muons) with momenta p > 0.3 (0.7) GeV/c in 
the laboratory frame. We merge photons consistent 
with e* decay radiation or bremsstrahlung with the 
corresponding e and veto events in which the e'^e~ 
pair is consistent with a photon conversion. We se- 
lect Kg candidates from tt+tt" final states for which 
the invariant mass is consistent with the nominal K'^ 
mass and the flight distance from the interaction point 
is larger than three times its uncertainty. We form tt*' 
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Figure 8: Predicted q^ dependence of the lepton forward- 
backward asymmetry in the SM (solid), the flipped-sign 
Cj model (dashed), the flipped-sign CgCio model (dotted) 
and the mirror image of the SM (dash-dotted) [ll| . 



candidates from two photons with energies larger than 
50MeV having an invariant mass of 115MeV/c^ < 






< 155MeV/c2. We require if* (892) candidates 
to have an invariant mass 0.82 < mK-rr < 0.97GeV/c^. 
Note that we imply charge conjugation throughout 
this article unless otherwise noted. 

We split the data set into two q^ regions, low q^ 
(0.1 < g2 < 7.02 GeVVc'') and high q^ (10.24 < q^ < 
12.96 GeVVc^ and g^ > 14.06 GeV'^/c'^). For the 
K*£~^i~ mode, we also report results that include the 
region of q^ < 0.1 GeV^/c^. We use two kinematic 
variables to select signal events, ttzes = a/s/4 — p'^ 
and AE = Eg - y/s/2, where p*g and E^^ are the B 
momentum and energy in the T(4S') center-of-mass 
(CM) frame, and ^/s is the total CM energy. We ex- 
tract signal yields from a one-dimensional fit to the 
TTiEs distribution for tties > 5.2GeV/c^, after a selec- 
tion on AE: -0.07 < A£; < 0.04 (-0.04 < AE^ < 
0.04) GeV for e+e^ [fji^ n^) events in the low q^ re- 
gion, and -0.08 < AE < 0.05 (-0.05 < AE < 0.05) 
GeV for e+e^ {fi^ n~) events in the high q^ region. 

The main backgrounds arise from combinations of 
leptons from two semileptonic decays {B —> Xit^v^ 
and B -> X2^-J>£ or S ^ Xl+ui and X -^ X't'Di). 
We suppress these combinatorial backgrounds by us- 
ing neural networks (NN). For each final state four 
separate NN are optimized to suppress either contin- 
uum or BB backgrounds in each of the two q^ regions. 
Inputs to these NN include event shape variables, ver- 
tcxing information and missing energy, where for each 
of the ten final states the NN selections are optimized 
to yield the highest statistical signal significance in the 
rriEs signal region {thes > 5.27GeV/c^). A potential 
background contribution arises from B -^ D{K'^*'^it)tt 
decays, where both pions are misidentified as muons. 
Therefore, we require the K^*>'k invariant mass to lie 
outside the 1.84—1.90 GeV/c^ region to veto this back- 
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ground. The final signal selection efficiencies vary 
from 3.5% for if+Tr^/i+^-at low values of q^ to 22% 
for K^ ■K~ e'^ e~ in the high q^ region. 

We use the vetoed J/ij] (7.02 < q^ < 
10.24 GeVVc"*) and ip{2S) events (12.96 < q^ < 
14.06 GeV^/c^) as control samples to validate our fit 
methodology, determine the probability density func- 
tions (pdf ) for the toes signal shapes and validate ef- 
ficiencies. Figure [S] shows the J /ip and ip{2S) branch- 
ing fractions we measure for the ten final states in 
comparison to the world averages [26l |. Our measure- 
ments agree well with the world averages. For all 
J/ip branching fractions and ipi^S) branching frac- 
tions measured in the e~^e~ mode, the uncertainties 
are smaller than those of the world averages. 
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Figure 9: Measurement of J/tp and ip{2S) branching frac- 
tions for the ten final states in comparison to world aver- 
ages [2g|. For the J/ip [^(25)] the test is performed in the 
low g^ and high q^ (entire) regions for fi'^ fJ,~ (solid points) 
and e"*'e~ (triangles) final states. 



We consider systematic uncertainties associated 
with reconstruction efficiencies, hadronic background 
parameterization in ^i^ IJ.~ final states, peaking back- 
ground contributions obtained from simulations, as 
well as possible isospin, CP, and lepton fiavor asym- 
metries in the background pdfs. We quantify the effi- 
ciency systematic using the vetoed J/ipK^*'> data sam- 
ples. These include charged track, 7r° and Kg recon- 
structions, FID, NN event selection, and AE and K* 
mass selection. The largest individual systematic er- 
rors come from hadronic FID, the characterization of 
the hadronic background and signal tties pdf shape. 
In the rate asymmetries most of these errors cancel at 
least partially. In general, the systematic uncertain- 
ties are very small compared to statistical uncertain- 
ties. 



3.1. Branching Fraction Measurements 

Figure [10] shows the TOes distributions for the 
K^,K'-\K*'^ and K*'^ modes after summing over 
e+e~ and iX^ IJ.~ modes, K*^ submodes and the low q^ 
and high q^ regions. Figure [Til shows the correspond- 
ing 771ES distributions in the low q^ region. We fit the 
TTiES distributions to extract signal and background 
yields, iVg and Nb respectively. We use an ARGUS 
shape '.21\ to describe the combinatorial background, 
allowing the shape parameter to float in the fits. For 
the signal, we use a fixed Gaussian shape unique to 
each final state, with mean and width determined 
from fits to the analogous final states in the vetoed 
J /tp events. We account for a small contribution from 
hadrons misidentified as muons by constructing a his- 
togram pdf using K^*'h n^ events weighted by the 
probability for the h^ to be misidentified as a muon. 
We also account for charmonium events that escape 
the veto, and for cross-feed contributions from misre- 
constructed signal events. In the entire q"^ region, we 
observe significant signal yields (> 4cr) in each of the- 
ses modes except for K''*i'^i~ where the significance 
is only l.Tcr. In the low q^ region significant yields 
(> 4cr) are only seen in the K+i+i' and K*+i+i- 
modes. The K'^£'^£~ mode is not seen at all. 
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Figure 10: Measured tties distributions in the entire q 
region (points with error bars) for K'^£^£~ (upper left), 
K°e+£- (upper right), K*+£+£- (lower left) and K*°e+£- 
(lower right) with fit results superimposed, full fit (blue 
solid curve), signal contribution (black dashed Gaus- 
sian), combinatorial background (magenta dashed curve), 
misidentified muons (green dotted histogram), and peak- 
ing backgrounds (red dotted Gaussian) . 



Averaging over i?+ and B^ modes in the low g^ 
region we measure preliminary branching fractions of 
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Figure 11: Measured tties distributions in the low q re- 
gion (points with error bars) for K'^t^l' (upper left), 
K°t+r (upper right), K'+l+T (lower left) and K*°e.+r 
(lower right) with fit results superimposed, full fit (blue 
solid curve), signal contribution (black dashed Gaus- 
sian), combinatorial background (magenta dashed curve), 
misidentified muons (green dotted histogram), and peak- 
ing backgrounds (red dotted Gaussian). 



results |2g|. The BaBarB{B -^ Ki+i') measurement 
lies m.ore than one standard deviation below the re- 
sults from Belle and CDF but agrees well with one SM 
prediction [ll|- For B ^ K* (.+£-, the BaBAR mea- 
surement is in good agreement with the CDF result 
and both SM predictions but lies more than one stan- 
dard deviation below the Belle result. 
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Figure 12: Preliminary BABAR B -^ Kt^l^ (open circles) 
and B ~» K*£^£~ (triangles) branching fraction measure- 
ments in the low q^ and high q^ regions. The dark-green 
shaded region shows the SM prediction for B -^ K*i^£~ . 



B(B^ Kt+r) = (1.8 ± 0.4 ± 0.08) X lO"'^ 
B{B^K*e+e') = (4.3tiJ±0.3) X 10"^ (6) 

that are of similar size as those in the high q^ region: 

B{B ^ K(+r) = (1.4 ± 0.4 ± 0.07) X 10"^ 
B{B ^ K*(+r) = (4.2 ± 1.0 ± 0.3) X 10"^. (7) 

Figure [12] shows the BaBAR measurements in com- 
parison to a SM prediction for K*i^i~ in the low q^ 
region [l^, [2g| . The isospin averaged branching frac- 
tion is in good agreement with this prediction. 

Combining the two q^ regions and correcting for the 
vetoed J/^ and '>p{2S) regions we measure preliminary 
total branching fractions of 



B{B^ Ki+r) 
B{B ^ K*e+£-) 



-7 



(3.9 ±0.7 ±0.2) X 10 
(ll.ltl:8±0.7) X lO^'^. (8) 



Figure [13] shows recent total branching fractions 
measurements for B -^ K£+£~ and B -^ K*£+£~ 
from BABAR [2i|, Belle [s^l, and CDF [HI in compari- 
son to two SM predictions [nl, [S^l- For completeness, 
we have also included previous B -^ Xst^i~ branch- 
ing fraction measurements from BaBAR 1331 and Belle 
[34| in comparison to a SM prediction |ll| . The new 
BaBAR exclusive measurements supersede the previous 
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Figure 13: Measured B -^ K^'''>£+£r and B -> Xs£+r 
branching fractions from BABAR (solid points). Belle (tri- 
angles) and CDF (squares) in comparison to the SM pre- 
dictions (shaded regions). 



3.2. Isospin Asymmetry Measurements 

Figure [H] shows the individual BaBAR B+ and 5° 
branching fraction measurements in the low q^ and 
high q^ regions as well as the SM predictions for the 
K*£+£- modes in the low q^ region [10, 28]. While in 
the high q^ region the isospin-related branching frac- 
tions for B+ -^ K'-*^+£+£- and B° -^ K'-*^°£+£- are 
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consistent, they differ considerably in the low q^ re- 
gion and deviate from the SM prediction by more than 
1(7. We see no signal events for B^ —^ Kgl'^l' and 
so we can only set a branching fraction upper limit 
at 90% confidence level. The preliminary BaBAR mea- 
surements in the low q^ region yield 



B{B° -> K"i+r) < 0.9 X 10"^ @90% C.L. 
B{B+ ^ K+e+r) = (2.5 ± 0.5 ± 0.1) X 10"'^ (9) 

B{B" ^ K*°£+£-) = (2.6ti:o±0.2) X lO"'^ 
B{B+ ^ K*+£+e-) = (9.8^2:4 ± 0.6) X 10-^.(10) 
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Figure 14: Comparison of B^ and B° branching frac- 
tion measurements for B —* K£^£~ and B — > K*£'^£~ 
modes in the low q^ and high q^ regions, where B{B^ -^ 
K+£+£-) [B{B° ^ K° £+£-)] are shown by light [dark] red 
triangles and B{B+ ^ K*+£+£') [B{B° ~* K*° £+£-)] by 
cyan open squares [blue open circles]. 

Table U and Figure [15] summarize our results for 
the isospin asymmetries that are obtained from direct 
fits to the K^*'>° and K^*^+ data samples taking into 
account the different 5° and B'^ lifetimes. Both in 
the high q^ and entire q^ regions, Aj is consistent 
with zero. In the low q^ region, however, we ob- 
serve significant isospin asymmetries. For Ke^e~ and 
K£^£~, the minimum lies in the unphysical region. 
For Kfj.^^^ due to large statistical uncertainties, the 
Af measurement in the low q^ region is consistent 
with both maximal isospin violation and isospin sym- 
metry. 

Defining statistical significance by N^^ = -\/2 A log £ 
where C is the likelihood function we determine at 
what values N„ the Aj — hypothesis is rejected 
by the data. Figure [16] shows the likelihood curves 
obtained from fits to the K£'^£^ and K*£+£~ data 
samples. The parabolic nature of the curves in the 
Ai > —1 region demonstrates Gaussian nature of our 
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Figure 15: Measured isospin asymmetries for B -^ K£'^£~ 
(triangles) and B -^ K* £'^£~ (open circles) in the low q^ , 
high q^ , and entire q^ regions. 



fit results in the physical region. The right-side axis of 
Figure \W\ shows purely statistical significances based 
on Gaussian coverage. Including systematic uncer- 
tainties, the significance of Ai in the low q^ region to 
differ from the null hypothesis is 3.2a for K£~^£~ and 
2.7a for K*£~^£'. We have verified these confidence 
intervals by performing fits to MC experiments that 
are generated with Ai = fixed. Using a frequentist 
approach we obtain results that are consistent with 
our above significance calculations. 

The highly negative Ai values for both K£^£~ 
and K*£^£r at low q'^ suggest that this asymme- 
try may be insensitive to the hadronic final state. 
Thus, we sum the K£+£'- and K*£+£- likelihood 
curves as shown in Figure [T^] from which we obtain 
Ai = -O.GAtoll ± 0.03 for the combined sample. In- 
cluding systematic errors we find a 3.9a significant 
deviation from the null hypothesis. If we Include 
the pole region {q^ < 0.1 GeV^/c"*) in K*e'^e~, the 
isospin asymmetry is reduced to — 0.25lg;^g ± 0.03. 
Given the large statistical error this result is con- 
sistent with the SM within two standard deviations. 
It is also interesting to compare the results with the 
isospin asymmetry measured in the K*"f modes. The 
BaBAR measurement of ^^ '' = 0.05 ± 0.058 agrees 
well with the SM prediction [3^ . The isospin asymme- 
tries determined for the charmonium control samples 
are plotted in Figure [18] and are in good agreement 
with the SM predictions. 

The isospin asymmetry for K*£'^£~ modes expected 
in the SM is shown in Figure [T7] For g^ -^ 0, the SM 
predict a positive value of 6 — 13%, which is opposite 
in sign to our observation. A model in which the sign 
of C7 is flipped provides qualitatively a better descrip- 
tion of our results than the SM. An SM calculation of 
K*~^ and K*'^ partial decay rates integrated over the 
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Table I Preliminary BABAR measurements of isospin 
asymmetries in different q^ regions. Uncertainties are sta- 
tistical and systematic, respectively. The last line shows 
K*e^e~ results including the q^ < 0.1 GeV^/c* region. 
Mode all q^ low q^ high q^ 



KtJ.+ti- 
Ke+e- 

Ke+e- 

K'e+e- 

K*e+e- 



,+0.29 
0.37 
f-0.39 



0.13lo;37±0.04 



-0.73+q;5p ± 0.04 
-0.37lo:3l ± 0.04 



-0.91 

-1.41^ 



+ 1.2 



±0.18 
±0.04 0.21"* 



a+0.35 



-1.43to;85 ±0-05 



-0.20 



-,+0.36 
-0.26 

+0.22 



-0.26 



+0.50 



±0.05 -0.08 



±0.03 



-0.12+°ll ±0.04 



-o.eetol? ± 0.02 

-0.56to.ll ±0.03 



_^^^ ±0.03 
0.28to;3o ± 0.03 
+0:2? ±0.05 



0.32to;l5 ± 0.03 



0.18"* 



±0.04 



K*e+e- 



-0.27^0 jg ± 0.03 



low q^ region yields an isospin asymmetry prediction 
of ^/ = -0.005 ±0.02 [i,[2|]. The BABmAf result 
is consistent with this prediction at the < 3a level. 
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Figure 18: Isospin asymmetry measurements for the char- 
monium control samples. 



Table II Preliminary Acp results. The uncertainties are 
statistical and systematic, respectively. 

Mode all g^ low 5^ high q'^ 






-CIS"* 



±0.01 -0.18 



0.02+jJ;2° ±0.02 



-0.23 



±0.02 0.17+°;24 ±0.02 
-0.18to;gg ±0.04 

K*e+e- o.oi1q;i5 ±0.01 o.oilo;2j±o.oi 0.09l°;2i ±0.02 



K*+e+e- 0.01 toi! ± 0.02 



+0.19 
-0.19 
+0.38 
0.38 
1.25 
0.24 



±0.01 



1+0.36 , 



Figure 16: Preliminary logarithmic likelihood curves for 
B -^ K£+£- (long-dashed curve) and B -^ K*£+r 
(short-dashed curve) modes using only statistical errors. 
The summed A log £ curves are shown by the solid curve. 



assume a common background ARGUS shape param- 
eter ^. The results for Acp are summarized in Table 
im We observe CP asymmetries that are consistent 
with zero as expected in the SM. 



T3 



< 

T3 




3.4. Tests of Lepton Flavor Asymmetries 
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q2[GeV2/c4] 



Figure 17: The q^ dependence of the isospin asymmetry 
for B -^ K*e+£~ in the SM and for the flipped-sign dr 
model [y|. 



Table Unl shows the BaBAR results for the lepton fla- 
vor ratios TZk and TZk' both including and excluding 
events with q^ < O.lGeV^/c^. The most significant 
deviation from the SM expectations is found in the 
low g2 region, where TZk = 0.40+°;^° ± 0.02 hes ~ 2a 
below the SM prediction of one. Figure [TOl shows TZk 
and TZk* measurements from BaBAR and Belle [3J] for 
the entire q^ region in comparison to the SM predic- 
tions [ll|. At the present level of precision, all results 
are consistent with lepton flavor universality. 



4. Angular Analysis of S ^ K*e'^t 



3.3. Direct CP Violation Measurements 

We extract the CP asymmetries by performing fits 
to the split B and B datasets in charge-conjugate fi- 
nal states using all modes except B — *■ Kg£~^£~. We 



The data selection for the angular analysis is per- 
formed in a similar way as that for the decay rates, 
except that some selection criteria are tightened to 
minimize backgrounds 36]. Using 384 million BB 
events we select the same six K*£^£r final states as in 
section [21 We measure Tl and Afb in the low q^ and 
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Table III Preliminary lepton flavor ratio results. The un- 
certainties are statistical and systematic, respectively. 

q^ Region 7?x* Rk 



all 



1.37lj;;|g ± 0.09 0.96tJ];34 ± 0.05 



all + Q^ < 0.1 1.101E;;32±0.07 
low l.Ol^ " 



low + g^ < 0.1 0.56; 



-0.44 
-0.29 
-0.23 
-1.42 
-0.78 



±0.04 



high 2.15lo 78 ± 0-15 l-Oetasi ± 0-06 



R. 



1=^ 



BaBar 349 fb'\ '08 

Belle 253 fb', '04 (preliminary) 

SM prediction 



Rk t^ 



1-^ 



R^iq'>0.1) 



0.5 



1 1.5 2 2.5 3 3.5 4 



Figure 19: Preliminary measurements of TZk and IZk' 
from BABAR (points) and Belle (triangles) in comparison 
to the SM predictions (black bars). 



the high q^ regions, where the low q^ region is reduced 
to 0.1 < g^ < 6.25 GeV^/c'* to ascertain that leakage 
from J /"ipK* background into the final data sample is 
rather small. The final reconstruction efficiencies for 
signal events vary from 1.5% for K^tt'^^^ ^~ in the 
low q^ region to 12.6% for K^n^e^e^ in the high q^ 
region. 

Due to the small event samples, we fit the data 
to the one-dimensional angular distributions in Eq[4l 
Note, that Tl is constrained by W^(cos 6'if) and for 
fixed J^L, Afb is constrained by W{cos9i). For each 
q^ region, we combine events from all six final states 
and perform three successive unbinned maximum like- 
lihood fits. First, we fit the mES distributions in the 
region tues > 5.2 GeV/c^ to obtain the number of 
signal (Ns) and background (Nb) yields using a Gaus- 
sian pdf for signal with mean and width determined 
from the vetoed charmonium sample and an Argus 
shape for the combinatorial background. We account 
for a small contribution from hadrons misidentified 
as muons, for misreconstructed signal events and for 
charmonium events that escape the veto. 

In the second fit to the cos 9k distribution for events 
with niEs > 5.27 GeV/c^ we extract !Fl, where 
the normalization for signal and background events is 
taken from the first fit to ruEs- In the third fit to the 
cos 9i distribution for events with niEs > 5.27 GeV/c^ 
we extract Afb for J-'l fixed to the result from the pre- 
vious fit and the normalization determined from the 
fit to niES- We model the cos 6'i<- and cos6'f shapes of 



the combinatorial background using e"*"e and fi'^ ij- 
events, as well as lepton fiavor violating e /x± events 
in the 5.20 < tues < 5.27 GeV/c^ sideband. The 
signal distribution is convolved with the detector ac- 
ceptance as a function of cos 9k and cos9i, respec- 
tively. The correlated leptons from B -^ D'^*^I!.v, 
D -^ K'^*^I!.v give rise to a peak in the combinatorial 
background at cos9i > 0.7 which varies as a function 
of ruEs- We consider this variation in our study of 
systematic errors. 

We test our fit methodology by using the large 
samples of vetoed charmonium events. Figure [^D] 
shows the cos 9k and cos 9g distributions for all B -^ 
J/i^K* events. We extract J^l = 0.569 ± 0.010 and 
Afb = -0.001 ±0.011. While Afb is consistent with 
zero as expected in the SM, J^l agrees with a recent 
BABAR measurement yielding J^l = 0.56 ± 0.01 [13] ■ 
The J^L and Afb results for the individual J/tpK* 
submodes are shown in Figures [211 and [22l We also fit 
the niEs and cos 6*^ distributions of the K^£^£~ de- 
cays and find Afb consistent with zero as expected. 
We have performed fits using signal events generated 
with different values of the Wilson coefficients cover- 
ing the allowed ranges of !Fl and Afb- We find no 
bias in extracting J-'j^ and Afb from the fits. 

Figure [23] shows the cos 9k and cos 9i distributions 
for the six K*£'^£~ modes combined in the low q^ and 
high q^ regions. The fits include pdf shapes of the sig- 
nal, combinatorial background, peaking backgrounds, 
and backgrounds from misidentified muons. In the 
low q^ and high q^ regions we measure 



jr|o» ^ 0.35 ±0.16 ±0.04 



r 



high 



and 



\ low 

y^FB 



A 



high 
FB 



0.7lto:22±0.05 



= 0.24l:°;^^ ± 0.06 
= 0.76^0.32 ± 0-07. 



(11) 



(12) 



Figures [M] show these results in comparison to the 
SM prediction and expectations of the three other 
models introduced in chapter [2] Though consistent 
with the SM prediction, the results for JF^ and Afb 
favor the flipped-sign Cj model [22] . Figure [25] shows 
our Afb results in comparison to the present Belle 
results [3g|. In the Belle analysis the low q^ region is 
divided into two bins where the high bin extends to 
8.7 GeV^/c'*. The high q^ region is divided into a bin 
between the J/4> and ip{2S) resonances and two bins 
above the ip{2S). The results of both experiments 
are in good agreement and the Belle results also favor 
the flipped-sign Cy model. The large values of Afb 
in the high g^-region disfavor the wrong-sign CgCio 
model at the > 3a level. Obviously, we need to per- 
form a model-independent global fit to measurements 
of Afb, ^l and other observables using data from 
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Figure 20; BABAR measurements of the cos 6k (top) and 
cos 6i (bottom) distributions (points with error bars) for 
the B -^ J/tp{-^ £'^£~)K* sample for all six modes com- 
bined with fit results superimposed, total fit (solid blue 
line), signal contribution (blue dots), combinatorial back- 
ground (green-dahed line) , self cross-feed (magenta dashed 
line) and feed-across background (red dashed line). 



BaBAR and Belle to look for significant discrepancies 
to the SM predictions and to extract non-SM contri- 
butions in the efi^ective Wilson coefficients Cr, Cg and 
CiQ. This, however, may require more precise results 
than are presently available. In addition, it would be 
useful to include measurements of other electroweak 
penguin decays such as B ^ Xs£^i^ in the fit. 
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Figure 21: BABAR measurements of J-l values for six 
J/tpK* control samples. The green bar represents the re- 
sult of a recently published BABAR analysis. 
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Figure 22: BABAR, measurements of Afb for six J/tpK* 
control samples. 



5. Search for B 



ni+t 



The exclusive decay B -^ tt£'^£~ is a. b ^ d tran- 
sition and is suppressed with respect to i? — > K£^£^ 
by |Vtd/Vtsp. The branching fraction expected in the 
SM is B{B+ -^ TT+£+£') ~ 3.3 x IQ-^ U,M- Belle 
has updated a search for B -^ n£^£^ for both tt^ and 
tt" in the recoil of an e'^e~ or /x"*"/.*" pair using 657 
million BB events ^0|- Since the main background 
originates from qq continuum events {q = u, d, s, c). 
Belle forms a Fisher discriminant from 16 shape vari- 
ables. They combine the Fisher discriminant together 
with the B vertex separation and the cos 9b distribu- 
tion into a likelihood ratio, where 6b is the angle of 



the B meson in the T(4S') rest frame with respect to 
the beam axis. After removing events in the vicinity 
of the J/tp and il}{2S) resonance regions they perform 
an unbinned maximum likelihood fit in the AE — tties 
plane. Since they observe no signal events, they de- 
termine branching fraction upper limits at 90% confi- 
dence level. The Belle limits are plotted in Figure [26] 
in comparison to BaBAR limits [41| that used 230 mil- 
lion BB events. Belle has set the lowest limits for 
TT+f+f" and for the combined 7r.^+^~ modes, whereas 
BaBAR has set the lowest limit for Tr'^£+£^. The Belle 
result of B{B+ -^ ir+f+fr) < 4.9 x IQ-^ @90% CL 
lies just a factor of ~ 1.5 above the SM prediction. 
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Figure 23; BABAR measurements of the cos 9k (top) and 
cos^f (bottom) distributions (points with error bars) for 
the six combined K*i^£~ signal modes in the low g^ (left) 
and high g^ (right) regions with fit results superimposed, 
total fit (red solid histogram), signal contribution (blue 
dashes), combinatorial background (black dots) and peak- 
ing backgrounds (green long dashes). 



6. Search for B -^ K^*^ 



vv 



The decays B — > K'^*^vv proceed throuf 
penguin and weak-box diagrams where U is 
pair. In the SM the branching fractions are 
to be [l|,[3|: 
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Figure 24: BABAR, measurements of Afb (top) and J-l 
(bottom) in the low (f and high q^ regions (points with 
error bars). For comparison, predictions are shown for 
the SM (solid curve), the flipped-sign C7 model (green 
dashed curve), the flipped-sign C9C10 model (magenta 
dotted curve) and the mirror image model of the SM (red 
dash-dotted curve). The solid (dashed) lines in the Tl 
plot show the SM (flipped C7 model) predictions after in- 
tegrating over the low q^ and high (f regions separately. 
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New physics may modify these predictions 
butions from new loop and box-diagrams 
particles in the loop (see Figure ^ may 
constructively or destructively with the sigr 
yielding enhanced or reduced branching frac 

BABMi has searched for B — > K*°iyD a 
K*+vv {B -^ K+vv) modes in the 
semileptonically tagged B -^ D'^*Hi' eve: 
454 (351) million BB events. The semilept 
are selected by constructing the angle 0b, Di 
the B and the D£ system from 
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Figure 25: Comparison of BABAR (open circles) and Belle 
(solid points) results for Afb as a function of q . 



COS0 



B,Dl — 



'^EbEbu 



•nir, — m 



Di 



2\pb\\pdi\ 



(14) 



where, E's, Eoi are the energies of the B and DI sys- 
tem, ruBjiTiDi are corresponding masses and pbtPdi 
are their momenta. For signal events, the cos9b,di 
distribution is bounded by the [—1,-1-1] interval, that 
may be slightly increased by resolution effects, while 



a large fraction of background events falls outside this 
region. The D'^ and D* mesons are reconstructed in 
several final states. In the recoil we select only K^ , 
K~^Tr~ , KgTT~^, and i^+7r° final states that have no ad- 
ditional tracks in the event. Due to the two missing 
neutrinos, we require a large missing energy. In the 
K~^iyD analysis signal selection is accomplished with a 
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multivariate analysis [42 
of merit \4M ■ 



that optimizes Punzi's figure 
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Figure 27: Measured Eextra spectra in B ^ K*uv modes: 
K*° -^ K+-K- (top left), K*+ -^ K%{-K+-K-)-K+ (top 
right), K*+ -^ K%{-k'^-k°)-k+ (bottom left), and K'+ -^ 
K'^-K (bottom right). Curves show the total fit (solid), 
signal (red dots) and background (hatched area). 



PUNZI 



^.,. 



na/2 



(15) 



where the significance is set to n^ — 3. In the K*vV 
analysis Punzi's figure of merit is used to optimize 
the selection criteria of six variables: cos 63,01, the 
ratio of second-to-zeroth Fox- Wolfram moments R2, 
the K* mass niK* , the lepton momentum p^ , the sum 
of the missing energy and missing momentum E^^^^ + 
'Vmiss^ ^'^'i the polar angle of the missing momentum, 
where the latter three variables are measured in the 
CM frame. We fit the extra neutral energy^ in the 
region 0.05 GeV < Eoxtra < 1-2 GeV to the expected 
signal and background shapes determined from signal 
and background MC samples, respectively. Figure [27] 
shows the observed Eextra distributions for the four 
K*vv modes. Semileptonic double tags (i.e. both B's 
decay to D'-*'>£i') are used in both analyses as a control 
sample to validate the simulation. 

We see no significant signal yield in any of these 
modes and set preliminary branching fraction upper 
limits at 90% confidence level of: 



= 4.2 X 10"^ @90% CL, 
= 9.0 X 10"^ @90% CL, 
= 21 X 10"^ @90% CL. (16) 



^energy of all neutral showers in the electromagnetic 
calorimeter not associated with signal or semileptonic tag. 
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Figure [551 shows the status of searches for exclusive 
{K,K* ,TT, p,(j))vi' modes from BaBAR. and Belle (for 
535 milfion BB events) [H]. The 90% CL branching 
fraction upper limits for Kvv {K*vv) modes are still 
a factor of four (seven) above the SM prediction. For 
B -^ K*vv the BaBAR upper limits are the lowest, 
while for _B+ -^ K^vv Belle has set the lowest upper 
limit of B{B -^ K+vv) < 1.4 x 10-^@90% CL. 
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Figure 28: Branching fraction upper limits at 90% confi- 
dence level for B —f hvv modes from BABMi (yellow bars) 
and Belle (dark-hatched bars), where the hadron {h) is 
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7. Conclusion and Outlook 

The B° and B+ averaged B -^ K^*H+i- total 
and partial branching fractions from BABAR agree well 
with the SM predictions. In the high q^ and entire 
q^ regions, isospin asymmetries are consistent with 
zero as expected in the SM. In the low q^ region, 
however, we observe large isospin asymmetries in the 
Ki+i~ and K*£'^i~ data. We find no evidence for 
Kg£^£~ events in the low q^ region. While for K£'^£~ 
the fit yields an unphysical central value we measure 
Ai = -0.56t oil ±0-03 for B -^ K*£+£-. This result 
is qualitatively more consistent with the flipped-sign 
Cj model than with the SM. Summing the Alog£ 
curves for K£~^£~ and K*£^£^ samples we measure 
an isospin asymmetry of Ai = —0.64^0 m ± 0.03 for 
the combined sample that excludes the null hypothesis 
with a significance of 3.9(7. If isospin asymmetry in the 
low q^ region expected in the fiipped-sign C-j model for 
the combined sample is similar or smaller than that 
for K*£^£^ , we may need additional sources of isospin 
violation to explain the BaBAR measurements. 

All observed CP asymmetries are consistent with 
zero. The measured Ti-xi*) ratios are consistent with 
lepton fiavor universality. The measurements of the 
K* polarization and the lepton forward-backward 
asymmetry are consistent with the SM prediction. 
Both jBaBAR and Belle results, however, favor the 
fiipped-sign Cy model over the SM. They disfavor 
the flipped-sign CgCio model at the > 3tT level. To 
proceed beyond these qualitative statements we need 
to perform a model-independent global fit to mea- 
surements of Afb, J'h and other observables using 
data from both experiments. In this fit we should 
include also corresponding measurements of the in- 
clusive B -^ Xs£'^£~ channel and results of other 
electroweak penguin decays, since these might help 
to improve precision in the extraction of the effec- 
tive Wilson coefficients. Presently, however, all exist- 
ing measurements and their SM predictions have too 
large uncertainties to uncover deviations from the SM. 
Eventually, a global model-independent fit may allow 
us to extract non-SM contributions (moduli and argu- 
ments) in the effective Wilson coefficients C7, Cg and 
Cio if they exist. 

Belle has set the lowest 90% CL branching fraction 
upper limit on _B — > vr+^+i?^ that lies about a factor 
of 1.5 above the SM prediction. BaBAR has set new 
branching fraction upper limits at 90% CL for B 



K+vD, B -^ K*°vD and B 



K 



*+, 



modes, BABAR. has set the lowest upper limits, while 
for B -^ K^vv the Belle upper limit is the lowest. 

The entire BABAR data sample consists of 465 mil- 
lion BB events. Though all present results in 
the exclusive modes will be updated with the full 
BaBAR data sample, the improvements in precision 
will be rather limited. We may succeed, however, in 



performing measurements of partial branching frac- 
tions, decay-rate asymmetries, longitudinal polariza- 
tion and lepton forward-backward asymmetry in four 
rather than two regions of q^. Furthermore, additional 
information will be obtained from an analysis of the 
sum-of-exclusive modes B — > KmT£~^£~ (n < 3) to ap- 
proximate the inclusive B — > Xs£~^£~ decays, in which 
we will measure the same observables with the full 
BABAR data set as in the exclusive analyses. A fully- 
inclusive B -^ Xst^£~ is presently not feasible, since 
the kinematic variables IS.E and tues are not defined 
here. In order to reject background from semileptonic 
decays an alternate strategy is needed. A powerful 
method consists of reconstructing one B meson com- 
pletely in a hadronic final state and then look for a 
lepton pair in the recoil. Since the efficiency for recon- 
structing hadronic final states is only ~ 10~^, a data 
sample much larger than presently available is needed 
to utilize this method successfully. 

Using the present analysis strategy we estimate 
about 15 billion BB events to measure ApB in 
K*£^£^ in eight bins of q^ with a precision of ~ 30% 
in each bin. For a similar analysis of the sum-of- 
exclusive decays about 7-8 billion BB events should 
be sufficient, if we reconstruct about half of all final 
states. A fully inclusive analysis, however, that looks 
for a dilepton pair from a common vertex in addition 
to fully reconstructed B in hadronic final states re- 
quires at least 50 billion BB events to see of the order 
of 100 events for an overall efficiency of 5 x 10~*. As- 
suming the SM branching fraction and an overall effi- 
ciency of 2.5 X 10"'', a sample of 2 billion BB events 



is needed to observe 10 events in B 



K' 



In 



similar size BB sample we expect about 10 tt'^£'^£~ 
events assuming the SM branching fraction and the 
performance of the previous BABAR analysis. 

The LHCb [Hi experiment at CERN will start data 
taking in 2008/2009 and the KEKB may be grad- 
ually upgraded to SuperKEKB [4y| in the coming 
years. LHCb and SuperKEKB will accumulate suffi- 
cient data to perform the Afb measurements in B —^ 
K*£^£~ and in the sum-of-exclusive B -^ Kmi£'^£~ 
modes. In addition they will observe B -^ tt£'^£~ . For 
the inclusive B -^ Xs£^£~ and B -^ K^*^iyp analyses, 
however, we need to wait for the high-luminosity Su- 
perB factory ^n\ that is proposed at Frascati. With 
a design luminosity of 10^^cm~^s~^ SuperB will pro- 
duce more than 10 billion BB events a year. 
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